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Fak, van Dijk and Wills (FDW) question our interpretation of elastic neutron-scattering exper- 
iments in the antiferromagnetic phase of UPt3. They state that our analysis is incorrect because 
we average over magnetic structures that are disallowed by symmetry. We disagree with FDW 
and reply to their criticism below. FDW alsopoint out that we have mistaken the magnetic field 
direction in the experiment reported in Ref. [i. We correct this error and note that our previous 
conclusion is also valid for the correct field orientation. 
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We disagree with the claim of Fak, et al.EI that our anal- 
ysis of elastic neutron-scattering experiments in the anti- 
ferromagnetic (AFM) phase of UPt3 is incorrect because 
we average over magnetic structures that belong to dif- 
ferent irreducible representations of the crystallographic 
space group. Classification of magnetic structures and 
magnetic phase transitions on the basis of irreducible 
representations of the space group and time-inversion ne- 
glects the fundamental role that pxjiange interactions 
play in magnetic phase transitionsBatrEl Exchange inter- 
actions are invariant under continuous rotations of all the 
moments, and typically dominate the anisotropy energies 
that couple the atomic moments to the lattice. Classifica- 
tion of magnetic structures based on the exchange group 
accounts for the wide variety of magnetic structures that 
are observed in magnetic materials. The Shubnikov clas- 
sification, which does not take into account the higher 
symmetry of the exchange interactions, disallows some 
of these structures.B 

Thus, for a magnetic instability driven by exchange 
interactions the primary irreducible representation is 
based on the combined group of continuous rotations in 
spin space, the crystallographic space-group and time- 
reversal, G ex . The irreducible representations of the ex- 
change group combine several irreducible representations 
of the space group.Q Thus, not only are magnetic struc- 
tures corresponding to irreducible representations of the 
space group allowed. On the contrary, structures that 
are a combination of irreducible representations of the 
space group, but belong to one exchange representation, 
are also possible magnetic structures. Many examples 
of magnetic structureS|-.with these type of "mixed spapft* 
group representations'^ are described in the literatureBEl 

In most materials the magnetically ordered phase is 
defined by one irreducible representation of the space 
group due to the anisotropy energies which resolve (at 
least partially) orientatkmal degeneracies within the ex- 
change representation.Ej However, since the anisotropy 
terms are relatively weak, the energy splitting of differ- 
ently oriented magnetic states are small. Thus, magnetic 
domain structures, including their response to magnetic 
fields, should be analyzed using the degenerate, or nearly 
degenerate, states within the full exchange multiplet. We 



believe this is the correct approach to understanding the 
magnetism and to analyze the possible magnetic struc- 
tures in the heavy fermion compound UPt3. 

In our analysis, we considered a general model for UPt3 
compatible with the available dataJlj We selected one ir- 
reducible representation of G cx that is consistent with 
elastic neutron scattering data in zero field. If we ne- 
glect the spin-lattice couplings then only the relative ori- 
entations of the atomic moments in the magnetic unit 
cell are fixed by the primary irreducible representation. 
Anisotropy energies are also included to resolve, or par- 
tially resolve, the degeneracies of the exchange represen- 
tation. 

Neutron scattering and X-ray experiments in UP^ 
show AFM order with propagation vector qi — aj/2.ll2l 
The magnetic U ions occupy two symmetry equivalent 
positions in the unit cell. The magnetic representation 
has 6 dimensions (3 times the number of magnetic ions). 
Until very recently, the crystal structure of UPt3 was 
thought to be hexagonal with space group D\ h . How- 
ever, a recent X-ray diffraction experiment revealed a 
lower trigonal symmetry with space group D\ d xs In ei- 
ther case, the magnetic representation can be decom- 
posed in six one-dimensional representations. Three of 
these correspond to FM alignment of the ions in the unit 
cell; the other three representations correspond to AFM 
alignments. The alignment of the magnetization or sub- 
lattice magnetization may be along the x, y or z axes. 
However, these six structures are connected with only 
two exchange representations corresponding to FM or 
AFM alignment in the unit cell. Table Q shows the ir- 
reducible representations and basis functions of the crys- 
tallographic space groups D\ h and D^ d grouped by their 
corresponding exchange multiplets. 

Our study is based on a free energy functional (Eq. 
9 of Ref. [ll]) which includes the exchange, anisotropy 
and Zeeman energies. First, a uniaxial anisotropy term 
(not shown in Eq. 9 of Ref. [ll]) restricts the order pa- 
rameter to the basal plane. In addition, the in-plane 
(hexagonal) anisotropy energy favors alignment of the 
moments along any of the three directions perpendicu- 
lar to the hexagonal lattice vectors. Note that the form 
of the anisotropy energy is the same for cither D^ h and 
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TABLE I: Irreducible representations and basis func- 
tions of the space groups D\ h and D\ d grouped by ex- 
change multiplets. FM and AFM refer to ferromagnetic 
or antiferromagnetic alignment of the two U ions on the 
unit cell. We use the notation of Kovalev in Ref. [L4|. 



that experiment the field was along the reciprocal lattice 
direction [-1,2,0]. The ratios reported in Eq. 5 of Ref. |Tl] , 
and in the paragraph that follows that equation, should 
be modified as follows. When only domain "1" is popu- 
lated we have r = 1. For a crystal with equally populated 
magnetic domains, the correct ratio between the scatter- 
ing rate at high field and zero field is 



1 - (0.441 cos((9 ff +tt/2)) 2 
(1 - (0.441 cos(6>)) 2 } 



= 0.89 . 



(1) 



D^ d symmetry groups. The effect of a magnetic field on 
the AFM order is included through the Zeeman coupling 
to the atomic moments, which in general mixes differ- 
ent nearly degenerate representations of the space group 
within the exchange multiplct.Eij The competition be- 
tween the anisotropy energy and the Zeeman coupling in- 
duces hexagonal modulations of the upper critical field as 
a function of the orientation of the field in the basal plane 
at the transition to the superconducting phaseJlSO The 
in-plane anisotropy energy is small, since a large in-plane 
anisotropy energy would produce an orthorhombic mod- 
ulation of the upper critical field, which is not observed. 
Higher order anisotropy termslifl which might resolve the 
remaining degeneracy and thus favor alignment of the 
moments along the propagation vector of the magnetic 
order would be extremely small. Therefore, the three 
structures shown in Fig. 1 of Ref. [ll] are degenerate, or 
quasi-degenerate, and certainly should be considered in 
the analysis of the magnetic structure and neutron scat- 
tering in the presence of an in-plane magnetic field. Thus, 
in our analysis we consider the possibility of degenerate, 
or nearly degenerate, magnetic structures by making an 
average over different distributions of domains. We also 
presented results and predictions for the single magnetic 
structure with the magnetization parallel to the propa- 
gation vector. The authors of the comment seem to have 
overlooked this prediction, which if we had confined our 
analysis to a single representation of the space group, as 
Fak, et al. advocate, would be the only relevant struc- 
ture. 

We did mistake the magnetic field direction in the ex- 
periment reported in Ref. 111. In the correct geometry of 



Our previous conclusion, stated for the incorrect field ori- 
entation, is unchanged for the correct field orientation; it 
is not possible based on existing data to conclude whether 
or not the U moments rotate with the field, because of 
the small change in intensity that is expected for this 
Bragg peak and the large error bars that are reported 
for the intensity. We also concluded that, in order to 
understand UPt3 magnetism in the presence of magnetic 
field or under pressure, systematic, zero-field measure- 
ments of the intensity of a number of magnetic peaks in 
the same single crystal, such as those reported in Ref. 
|l9| , need to be carried out. Furthermore, our hypothesis 
that intrinsic stacking faults pin the AFM domain walls 
in the ab-plane and fix the spatial distribution of domains 
with different propagation vectors has been recently rein- 
forced. For uniaxial pressures applied to the basal plane 
a significant increase in the magnetic intensity has been 
reportedEj in contrast with the relatively small change in 
a magnetic field.BEil Pinning by intrinsic stacking faults 
may help explain this difference, since the applied mag- 
netic field leaves the distribution of regions with different 
propagation vectors unaltered. However, uniaxial pres- 
sure likely disturbs the configuration of stacking faults 
leading to a stronger effect on the magnetic structure. 

In conclusion, our analysis of the neutron scattering 
data is based on a sound theoretical model for possi- 
ble magnetic structures in UPt3, which is more general 
than would be allowed based on a single irreducible rep- 
resentation of the space group. The relative importance 
of exchange interactions leads naturally to mixed irre- 
ducible representations of the crystal space group, which 
are relevant because they are energetically allowed. 



* Present address: Department of Physics and Astronomy, 
Clemson University, Clemson SC 29634 

1 N. H. van Dijk, B. Fak, L.P. Regnault, A. Huxley and 
M-T. Fernandez-Di'az, Phys. Rev. B 58, 3186 (1998). 

2 B. Fak, N. H. van Dijk, and A. S. Wills, Phys. Rev. B 
(2002), to appear. 

3 J.-C. Toledano and P. Toledano, The Landau Theory of 
Phase Transitions (World Scientific, Singapore, 1987). 

4 I. E. Dzyaloshinskii, Sov. Phys. JETP 19, 960 (1964). 

5 A. F. Andreev and V. I. Marchenko, Sov. Phys. JETP 43, 
794 (1976). 



6 A. F. Andreev and V. I. Marchenko, Sov. Phys. Usp. 23, 
21 (1980). 

7 Yu. A. Izyumov, V. E. Naish and R. P. Ozerov, Neu- 
tron diffraction of magnetic materials (Consultants Bu- 
reau, New York, 1991). 

8 For some materials, such as the orthoferrites, approxi- 
mately 20% of the magnetic phases are described bm mul- 
tiple irreducible representations of the space group.u 

9 A. Oles, F. Kajzar, M. Kucab and W. Sikora, Magnetic 
structures determined by neutron diffraction (Warszawa, 
Krakow, 1976). 



3 



10 Anisotropy energies, which arise from spin-orbit, dipolar 
and indirect interactions, may also lead to mixing of mag- 
netic structures belonging to different irreducible represen- 
tations of G ex . 

11 J. Moreno and J. A. Sauls, Phys. Rev. B 63, 24419 (2001). 

12 The propagation vectors §2 = a%/2 and cfa = (a* — a|)/2 
are also present since the wave vector has three arms. 

13 D. A. Walko, J.-I. Hong, T. V. C. Rao, Z. Wawrzak, D. N. 
Seidman, W. P. Halperin, and M. J. Bedzyk, Phys. Rev. 
B 63, 054522 (2001). 

14 O. V. Kovalev, Irreducible Representations of the Space 
Groups (Gordon and Breach, New York, 1965). 
Anisotropy energies can also couple an AFM representa- 
tion with the FM order parameter belonging to a differ- 
ent exchange multiplet. This leads to weak, Dzyloshinski- 
Moriya ferromagnetism. Since the induced magnetization 



couples linearly to the external field we included this cor- 
rection fn-aur analysis, in addition to the quadratic Zeeman 
coupling.EJ. 

16 N. Keller, J. L. Tholence, A. Huxley, and J. Flouquet, 
Phys. Rev. Lett. 73, 2364 (1994). 

17 J. A. Sauls, Phys. Rev. B 53, 8543 (1996). 

18 For example, terms resulting from the coupling of different 
exchange representations. 

19 A. I. Goldman, G. Shirane, G. Aeppli, B. Batlogg, and 
E. Bucher, Phys. Rev. B 34, 6564 (1986). 

20 N. H. van Dijk, P. Rodiere, F. Yakhou, M.-T. Fernandez- 
Dfaz, B. Fak, A. Huxley, and J. Flouquet, Phys. Rev. B 
63, 104424 (2001). 

21 B. Lussier, L. Taillefer, W. J. L. Buyers, T. E. Mason and 
T. Petersen, Phys. Rev. B 54, R6873 (1996). 



